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The absence of a light Higgs boson causes vector boson couplings to become strong at 1 TeV. 
A general framework for a systematic and consistent treatment is provided by effective theories 
of electroweak symmetry breaking. Already in next-to-leading order there appear quartic gauge 
couplings that go beyond the standard model and are hence called anomalous. We investigate 
intermediate three gauge boson states W + W~Z and ZZZ occurring in six quark production in 
electron positron collisions under the conditions of the International Linear Collider. We perform 
a sensitivity analysis of the relevant anomalous quartic gauge couplings presenting their expected 
limits. 



PACS numbers: ll.15.Ex, 11.30.Qc,12.15.-y, 12.39.Fe 

I. INTRODUCTION 

One important objective of future accelerator exper- 
iments is to unveil the mechanism behind electroweak 
symmetry breaking (EWSB). In this respect generic ef- 
fective field theories are a proper choice for a systematic 
and consistent treatment. These theories arise in a natu- 
ral way, if no light Higgs is assumed to be present in the 
particle spectrum Including all operators obeying 
the underlying symmetries, these effective theories are 
non-renormalizable in a narrow sense. At each order of 
perturbation theory more counterterms have to be intro- 
duced in order to render the whole theory finite. These 
so called anomalous couplings parameterize the effects of 
new physics in a generic, model independent way. Since 
the construction of the theory starts from the known low 
energy behavior of the electroweak interaction, the way 
of describing EWSB is called the bottom-up approach, 
see e.g. ref. [2J. Alternatively, the top-down approach 
constitutes the whole particle spectrum and symmetries 
of the theory at high energies from the beginning by con- 
structing renormalizable Lagrangians. 

In the following we consider the interaction between 
four bosons, in particular production of intermediate 
three boson states. The relevant CP conserving opera- 
tors arising in next-to-leading-order are the two SU(2)y- 
conserving Q 
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gauge-covariant derivative. The field £ parameterizes the 
Goldstone sector and provides a generally nonlinear real- 
ization of the electroweak gauge symmetry 0] . In unitary 
gauge £ = 1 the operator reduces to 
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where and r 3 are the generators of weak isospin group 
SU(2). Three further operators contributing to a generic 
four-point-interaction are 
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which, in contrast to Eqs. (UKU), violate 5/7(2) via the 
appearance of the operator T = Xr 3 X^. 

These interaction terms are evaluated in three-boson- 
production via the processes 



(W + W~Z) 
(ZZZ) 



> qqqqqq 
qqqqqq. 
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Only the sum of 5t/(2)-conserving and -violating op- 
erators contributes to each of these processes. Hence, 
there will be no direct sensitivity to S l L7(2)-violation. 
This question can be accessed by complementary anal- 
ysis, such as W-scattering, see e. g. 0,[j| . The coupling 
parameters can be mapped to masses of resonant par- 
ticles that may be indirectly accessible in the range of 
future colliders 0, Q . 

The results presented here are based on previous anal- 
yses In this paper we now extend this work in two 
directions: First, we consider the influence of the polar- 
ization states of the bosons. This is particularly interest- 
ing since the longitudinal polarized components directly 
correspond to the would-be Goldstone bosons associated 
with the symmetry breaking pattern. Second, the re- 
striction of considering on-shell bosons in the final state 
is abandoned in favor of investigating the complete set of 
six quark final states, corresponding to the pure hadronic 
decay of the intermediate bosons, which brings the over- 
all analysis closer to the future experimental situation. 

For a motivation of how the helicity angle may improve 
sensitivity we show Fig. [H presenting the reconstructed 
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FIG. 1: Reconstructed differential cross sections normalized 
to the SM cross section. 
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helicity distribution in dependence of the anomalous cou- 
pling parameters. We have chosen ct4 — o.§ — 0.1 and 
compare the helicity distributions implied by taking into 
account the anomalous couplings to the one of the stan- 
dard model. The distribution for cosQ* z are given after 
reconstruction as explained in the next sections. 

To guide the eye the histogram has been fitted by 
a simple second order polynomial. From inspection of 
Fig. [1] it is obvious that the appearance of anomalous 
couplings change the shape of the distribution from a 
rather flat one to a more pronounced distribution. Al- 
though the helicity state of the intermediate boson is not 
reconstructed this way it is possible to take into account 
the effects of the intermediate state of polarization. 

The helicity angle is given with respect to the rest sys- 
tem of the decaying boson, in this case the Z boson in 
WW Z intermediate states. To this end the co-ordinate 
system is properly boosted utilizing the reconstructed 
momenta (subject to kinematical fit). 

The curves in Fig. [1] are normalized to the respective 
total number of events to separate effects on total cross 
section from variations in the distribution itself. 



II. ANALYSIS LAYOUT 

We analyze intermediate three-boson-production at a 
center-of-mass energy of y/s = 1 TeV. The total inte- 
grated luminosity is assumed to be J C — 1000 fb _1 . 
An earlier analysis has shown Q that the sensitivity in- 
creases remarkably by using polarized initial e + e~-states. 
We therefore focus on the situation with 80 % right 
handed polarization of electrons and 60 % left handed 
polarization of positrons that has been the most favor- 
able scenario. The degree of polarization has been carried 
over from our earlier analysis. 

We consider all the six quark final states listed in 
Tab.H 

Each of the flavor final states is equivalent by a large 
number of Feynman graphs. They contain the signal 
events as well as a huge amount of background arising 



TABLE I: Generated six quark production processes at y/s — 
1 TeV considered in the analysis. (Cross sections including 
numerical errors are given). 
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from (a few thousand) standard model graphs with the 
same in and out states but without any quartic gauge 
couplings. (Note, however, that the actual generator al- 
gorithm differs from evaluating Feynman graphs, but us- 
ing technique described closer in [3]) 

The dominating part of this background stems from 
direct top production, via 

e + e~ — > tt > qqqqbb, (8) 
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which manifests itself as a six-jet-topology. They are 
indicated as tt in Tab. |TJ 

An additional type of background is related to misiden- 
tification. In case of the ZZZ intermediate states, also 
the contamination of reconstructed Zs with misidentified 
Ws has to be taken into account. Particle misidentifica- 
tion amounts to 8% of all reconstructed Zs. The converse 
case (Z identified as W in reconstructed WWZ events) 
is less important, since the ZZZ production is largely 
suppressed compared to WWZ production. 

The generation of events is achieved by using 
Whizard and Pythia @ is utilized for fragmentation 
to hadronic final states. Initial state radiation (ISR) has 
been taken into account. In order to reduce statistical 
fluctuations, the number of analyzed events is increased 
by a factor of hundred of the expected number (assuming 
a luminosity of C = 1000 fb _1 at the ILC) and finally 
renormalized. This procedure leads to a more reliable 
figure of merit for the sensitivity. Subsequent detector 
simulation is done with the fast simulation tool Simdet 
0- 



A. Fit Method 

We use a x 2_n t to extract information on anomalous 
couplings from the distribution of kinematical variables 
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The indices ikimn stand for the discretized phase space 
variables, which are six in number, and are used to de- 
scribe the event. They will be dropped in the following 
formulas for convenience. Confidence levels for a 4 and a§ 
are calculated with Minuit [ll[. The kinematical vari- 
ables, the distributions of which are sensitive observables 
to anomalous quartic gauge couplings and are therefore 
used in the analysis, are 



• the invariant mass My/z 
W\ - Z-subsystem 



V(pwi +Pz) 2 of the 



• the invariant mass My/w — \J {pwi + Pw 2 ) 2 of the 
Wi-W2-subsystem 

• the cosine cos# of Z momentum w.r.t. the beam 
axis 

• the helicity angle of each boson: 9yy ,#yy 2 , d* z . 

Since the bosons are off-shell, there are no self-energy 
corrections due to AQGC at NLO. Such contributions 
only arise at NNLO via loop-graphs including at least one 
AQGC- vertex. So in the considered order of perturbation 
theory there are no changes in the energy-momentum- 
relations compared to the SM, hence off shell effects do 
not lead to additional kinematical variables. 



The cross section der NLO in NLO is of second order in 
the coupling constants since the processes ([7]) contain 
at most one vertex with AQGC, 

da NLO = da SM + da^Oi + ^ daf e a k a e . (10) 

i kl 

Here d<r SM denotes the standard model contributions, 
der lnt the interference term linear in a, and dcr ac the 
quadratic term containing only AQG-vertices. The full 
information on the dependence on anomalous gauge cou- 
plings is incorporated in at most 5 different event depen- 
dent parameters, viz. 
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where the parameters A, . . . , E depend on the kinemat- 
ical variables, but not on anomalous couplings. They 
can be obtained by re- weighting the SM matrix elements 
at five different points in the a^-as-plane for a fixed 
event sample and solving the respective system of lin- 
ear equations obtained from Eq. (fTTjl . We choose 014 = 
±0.1, 0:5 = 0; a 4 = 0,0:5 = ±0.1, and a 4 — 0:5 = 0.1. 
The resulting \ 2 is an analytic function of oti due to (jTTJ) 
and can easily be minimized. 



B. Particle Identification and Cuts 

All processes used in the analysis are listed in Tab. [JJ 
In the following, for convenience and also historic rea- 
sons, processes that contain intermediate tt-production 
are referred to as "background" , whereas those contain- 
ing intermediate boson states only are termed "signal" 
events. In fact kt production vastly dominates the triple 
boson production by two orders of magnitude (see Tab. 
[TJ) . Note that a negligible fraction of three boson pro- 
duction is also included in the "background" , since they 
cannot (at generator level) be separated from each other 
in a six quark production. 

Due to hadronization the six quark final states form 
jets. Jet identification is done using the Durham al- 
gorithm [To| . In addition, we require a cut on the jet 
energy £j ct > 5 GeV and an upper bound on missing 
energy and momentum E^ is + p\ < (65 GeV) 2 to ac- 
cept only well defined jets. The reconstructed jet mo- 
menta are pairwise combined to form invariant masses 
-Mjetpairs = a/ (pji + Pj2) 2 of the boson candidates. Tops 
are reconstructed via b-decay combining 3 jets and re- 
quiring two of them to form a W . Cuts on respective 
invariant masses used to identify vector bosons and top 
quarks are given in Tab. [Til an d Fig. [2] Invariant masses 
of accepted candidates are required to have a minimum 
deviation from the nominal masses My/ = 80.403 GeV, 
M z = 91.1876 GeV, and M t = 174.2 GeV as given in 



For background identification we cut on the energy E 
and momentum p distribution of the top candidates in 
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particle 



lower limit upper limit 



W M w - 7.5 GeV (M w + M z )/2 
Z {Mw + M z )/2 M z + 7.5 GeV 
t M t - 15 GeV M t + 15 GeV 

TABLE II: Mass cuts used for particle reconstruction. 



10000 





FIG. 2: The upper figure shows the invariant mass of jet 
pairs from which the bosons are reconstructed for back- 
ground (red) and events without tt-production (black). The 
vertical lines indicate the selection parameters Mw = 
72.9 GeV . . . 85.8 GeV, M z = 85.8 GeV . . . 98.7 GeV. Tops 
are reconstructed by combining the best W in the event with 
all remaining jets and selected as indicated by the vertical 
lines in the lower figure \M t - M? om \ < 15 GeV. 
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(a) Energy distribution of signal (black) and background 
(red). 
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(b)Momentum distribution of signal (black) and background 
(red). 

FIG. 3: The cuts applied to top candidates after detector sim- 
ulation are shown: |£-486 GeV| < 55 GeV, |p-452 GeV| < 
63 GeV. These numbers are mean values obtained from 
Monte Carlo simulation of background events e + e~ — > ti. 
The plots show (a) the E and (b) the p distribution from 
tt background (red curve) and from the WW Z signal (black 
curve) . 
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the event as shown in Fig. [3] As "signal" events, only 
those are included that do not contain any contributions 
from tt-production, whereas the "background" consists 
to the largest part of such events. 

To enforce conservation laws a kinematical fit is ap- 
plied, which minimizes the missing energy and momen- 
tum constrained by the reconstructed mass, i.e. the fol- 
lowing expression for the reconstructed particles (reco) is 
minimized 
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III. SENSITIVITIES 

As a measure of sensitivity we use the 1-er standard de- 
viations of the a parameters obtained from minimizing 
X 2 of Eq. (J9j) with a SM distribution as null hypothesis. 
Here, we implicitly assume that the standard deviations 
(sensitivities) do not depend on the actual values of the 
anomalous couplings. This is reasonable as long as the 
true values for the anomalous couplings are unknown and 
are at least expected to lie within the range where the ef- 
fective theory approach is feasible. The resulting contour 
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FIG. 4: Obtained sensitivity for reconstructed events. Inner 
and outer contour delineate 68 % and 90 % confidence level, 
respectively. 

gives an estimation of the maximum probability that a 
measured event distribution is compatible with the null 
hypothesis. The sensitivities obtained by this procedure 
are shown in Fig. 3] for intermediate WWZ and inter- 
mediate ZZZ events, respectively. The combined \ 2 fit 
using reconstructed WWZ and ZZZ events can be found 
in Fig. [5l Tab. IIIII shows the 1-a standard deviation for 
the coupling parameters calculated from Minuit. 
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TABLE III: 1-a standard deviations for 04 and as obtained 
from reconstructed WWZ, ZZZ events and from a combined 
fit. The lower table shows the previously obtained results 
(taken from 
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FIG. 5: Combined fit from reconstructed WWZ and ZZZ 
events. 
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IV. CONCLUSION 

Based on our previous analysis the sensitivity on quar- 
tic gauge couplings of processes at the ILC involving 
three intermediate bosons have been improved. The im- 
provement has been in two directions: First we have in- 
cluded observables related to the polarization of the in- 
termediate bosons. This leads to a better sensitivity com- 
pared to the previous results. Second, we have utilized 
full six quark final states in the event generation that be- 
came available only recently. Matrix elements of all six 
fermion final states that can be accessed in the 1 TeV 
range have been included. This leads to a more real- 
istic description of the processes closer to the expected 
experimental behaviour. It also leads to a more realistic 
treatment of background events compared to former anal- 
ysis where only non-interfering tt-production has been 
assumed. 



Of all possible generic quartic gauge couplings the pro- 
duction of intermediate three boson states WWZ and 
ZZZ is sensitive to the ones involving a^fa%, and a^/aj. 
The present analysis is done for the ILC with the option 
for initial state polarization. Due to the large number 
of events that have to be considered, we have merely fa- 
vored a scenario of fully polarized initial states that has 
been the most sensitive case in our previous analysis [f| 

Presently, we have only utilized hadronic decays. 
These amount to 32% of the total events with intermedi- 
ate three bosons. A more complete investigation includ- 
ing semileptonic decays, i.e. e+e~ — * (VVZ) — > Aq + ££ 
is left for future work. 

A next task is to combine the present analysis with 
the results from IFW-scattering along the lines of [j| and 
update the characteristics of new resonant states that can 
be inferred from the present work. 
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